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Monte Carlo simulations are used to study the translocation of a polymer into and out of a ellipsoidal cavity 
through a narrow pore. We measure the polymer free energy F as a function of a translocation coordinate, 
s, defined to be the number of bonds that have entered the cavity. To study polymer insertion, we consider 
the case of a driving force acting on monomers inside the pore, as well as monomer attraction to the cavity 
wall. We examine the changes to F{s) upon variation in the shape anisometry and volume of the cavity, the 
polymer length, and the strength of the interactions driving the insertion. For athermal systems, the free 
energy functions are analyzed using a scaling approach, where we treat the confined portion of the polymer to 
be in the semi-dilute regime. The free energy functions are used with the Fokker-Planck equation to measure 
mean translocation times, as well as translocation time distributions. We find that both polymer ejection and 
insertion is faster for ellipsoidal cavities than for spherical cavities. The results are in qualitative agreement 
with those of a Langevin dynamics study in the case of ejection but not for insertion. The discrepancy is 
likely due to out-of-equilibrium conformational behaviour that is not accounted for in the FP approach. 


I. INTRODUCTION 

Polymer translocation is a fundamental process in 
which a polymer is transported through a narrow hole 
in a barrier divides two separate spaces Over the past 
two decades, there has been considerable progress in 
the development of techniques for detecting and mon¬ 
itoring single-molecule translocation events. Much of 
this work has been motivated by the promise of an ef¬ 
ficient and accurate translocation-based method for nu¬ 
cleotide sequencingOther technological applications 
include protein analysis,— filtration of macromolecules,— 
and controlled drug delivery!^ Polymer translocation is 
also an important part of numerous biological processes, 
including viral DNA packaging and ejection, transport 
of mRNA through the nuclear pore complex, horizon¬ 
tal gene transfer between bacteria, and protein transport 
across biomembranes— Due to its wide range of appli¬ 
cations, polymer translocation has been the subject of 
numerous theoretical and computer simulation studies in 
recent years. Much of this work has been summarized in 
several recent reviews 

One important type of polymer translocation involves 
movement of polymers into or out of confined spaces. 
Recent theoretical and computer simulation studies in 
this area have mainly focused on confinement in spher¬ 
ical or ellipsoidal cavitiesi^^— or laterally unbounded 
spaces between flat walls Much of this work is mo¬ 
tivated by the problems of viral DNA packaging and 
ejection, in which DNA is confined to a space with di¬ 
mensions comparable to that of its persistence length 
at near crystalline densities and very high internal pres¬ 
sures. Experimental studies suggest that these processes 
do not follow simple quasistatic dynamics. For example, 
ejection proceeds in rapid transient bursts separated by 
pauses^ while ultra-slow relaxation and nonequilibrium 
dynamics has been observed during packaging—^ A de¬ 
tailed understanding of these processes using theoretical 


methods requires, at a minimum, the use of semi-flexible 
chain models to account for the high energetic cost of 
confinement. Numerous studies have examined such 

On the other hand, such 
studies are complemented by those that use flexible-chain 
models to elucidate the specific role of conformational en¬ 
tropy on translocation.— other 

studies have examined the effects of solvent quality^, 
electrostatic interactions,—!^, temperature,— and ad¬ 
sorption to the cavity surface.— >21 

A few theoretical studies have considered the effect 
of cavity shape anisometry on polymer insertion and 
ejection— Ali et al. compared packaging and ejec¬ 
tion in ellipsoidal and spherical cavities of equal volume 
and found that flexible polymers package more quickly in 
spherical cavities but eject faster in ellipsoidal cavities— 
By contrast, both processes are faster for spherical cav¬ 
ities in the case of semiflexible polymer, leading the au¬ 
thors to suggest this as a reason for the spherical shapes 
of viruses with pressure-driven ejection. Zhang and Luo 
recently used a 2-D system to study translocation of a 
polymer into an elliptical cavity.^^L^ For flexible chains, 
they found that the translocation time increased with 
cavity anisometry— In the case of semiflexible polymers, 
they found that the packaging rate depended the location 
of the entry point to the cavity, as well as the packing 
fraction and chain stiffness— For a given elliptical cav¬ 
ity at high confinement, they found that entry along the 
semi-minor axis gave the fastest packaging for sufficiently 
high stiffness.— 

In many simulation studies of polymer translocation 
into or out of confined spaces, the scaling of the rate of 
transport with polymer length, N, and cavity dimension, 
R have been measured. Typically, the results are inter¬ 
preted using estimates of the confinement free energy, 
F, of the polymer in the cavity and how it varies with 
the degree of translocation. The case of polymer ejection 
from spherical cavities provides a noteworthy illustration. 
In an early study on the topic, Muthukumar assumed a 
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scaling F ^ N/R^/'' to describe the confined portion of 
the chain, where v k, 0.588 is the Flory scaling exponent. 
This was used to explain the observed scaling of the mean 
exit time of r , where (/> is the packing 

fraction of the cavity for a fully inserted polymerSub¬ 
sequently, Cacciuto and Luijten noted that the appro¬ 
priate scaling of the free energy for triaxial confinement 
is F ~ . Using an approach suggested by 

the scaling of r for driven translocation given in Ref. 
they predicted r ^ which was consistent 

with Monte Carlo (MC) dynamics simulationsSakaue 
later questioned the implicit assumption of a constant 
free energy gradient during translocation and derived an 
alternative scaling relation that accounts for the decrease 
and eventual loss of entropic driving force near the end of 
the processi^ The scalings observed in each of the MC 
dynamics studies of Refs. M and [l^ were shown to be 
limiting cases of this more general result. 

An alternative approach to using analytical estimates 
of the translocation free energy function is to calculate 
it explicitly for a chosen model using simulations. Ras¬ 
mussen et al. carried out such calculations for translo¬ 
cation of Lennard-Jones chains into absorbing spherical 
cavities used the Incremental Gauge Cell MC method.— 
The free energy functions were then used in conjunc¬ 
tion with the Fokker-Planck (FP) equation to calculate 
translocation time distributions and probabilities. They 
observed an interesting non-monotonic dependence of 
translocation times with a sharp peak located at a lo¬ 
cal free energy minimum.^ These results were consistent 
with other results by the same group obtained using self- 
consistent field theory.— The validity of the FP approach 
requires quasistatic conditions during translocation. As 
noted by Kantor and Kardar— this condition cannot be 
satisfied for long polymer chains. In addition, nonequi¬ 
librium behaviour has been observed in simulations of 
polymer ejection— On the other hand, for sufficiently 
high pore friction quasistatic translocation is possible 
and the FP approach is valid.— Recently, we used a MC 
method to calculate translocation free energy functions^^ 
and showed that the predicted translocation time dis¬ 
tributions were perfectly consistent with those obtained 
from dynamics simulations when the pore friction was 
sufficiently high— 

The present work is a theoretical study of the translo¬ 
cation of a flexible polymer into and out of an ellipsoidal 
cavity. Following our other recent workj ^^i^^d^ we use 
MC simulations to measure the translocation free energy 
functions. In addition to entropic effects, we also con¬ 
sider the effects of a force located in the nanopore that 
drives polymer insertion into the cavity. We also con¬ 
sider the effect of monomer attraction to the cavity wall, 
as in Ref. [13 We examine the effects on the free energy 
of varying several key system parameters, with special 
attention to the effects of the anisometry of the cavity. 
The free energy functions are used with the FP formalism 
to calculate translocation times for both polymer inser¬ 
tion and ejection. We find that ejection is predicted to be 


faster in ellipsoidal cavities in agreement with Ref. [l^. In¬ 
terestingly, we also find that insertion tends to be slower 
for ellipsoidal cavities, in disagreement with results from 
that study. This points to the importance of nonequi¬ 
librium dynamics, which is not accounted for in the FP 
approach, as well as the limitations of using free energy 
functions to predict translocation dynamics. 


II. MODEL 

We employ a minimal model of a polymer chain that 
translocates through a narrow pore between an ellipsoidal 
cavity and semi-infinite space on one side of a flat wall. 
The polymer is modeled as a flexible chain of N hard 
spheres, each with a diameter of cr. The pair potential 
for non-bonded monomers is thus Mnb(?') = for r < 
tj, and Unb(?’) = 0 for r > cr, where r is the distance 
between the centers of the monomers. Pairs of bonded 
monomers interact with a potential Uhir) = 0 if 0.9cr < 
r < 1.1(7, and u\^{r) = oo, otherwise. Consequently, the 
bond length can fluctuate slightly about its average value. 
Each monomer interacts with the walls of the system 
(pore, cavity and barrier) with a hard-particle potential. 
Thus, the monomer-wall interaction potential is Mw(^) = 
0 if the distance r between the center of the monomer 
and the nearest point on the wall satisfies r > 0.5cr, and 
Mwaii = oo if r < 0.5cr. In addition, we consider two 
different potentials designed to drive the polymer into 
the cavity. In the first case, monomers that lie inside the 
pore are subject to a potential described by a constant 
driving force of magnitude /d that is directed toward the 
cavity. In the second case, monomers inside the cavity 
whose centers lie a distance less than cr from the wall 
have a potential energy of —e (where e > 0), while all 
other monomers have zero potential energy. 

The pore connecting the cavity to the open space is 
cylindrical in shape with a diameter D and length L. 
For most simulations in this study, we choose D = 1.2a 
and L = l.Scr. We choose the z axis to pass through the 
center of the cylindrical pore. The semi-infinite space 
on one side of the pore is bounded by an infinite flat 
wall perpendicular to z. The cavity is an ellipsoid of 
revolution, with a semi-axis length of a along z and b 
along X and y. We consider the cases of both prolate 
(a > b) and oblate (a < b) ellipsoids, as well as the 
special case of a spherical cavity {a = b). The volume V 
of the ellipsoid is defined to be the volume of the subspace 
accessible to the centers of the monomers. This subspace 
is enclosed by a virtual surface, each point on which lies a 
distance of 0.5cr from the nearest point on the cavity wall. 
This virtual surface deviates somewhat from ellipsoidal 
geometry, except in the special case where a = b. The 
aspect ratio, r, of the cavity is defined to be the ratio of 
the dimensions of the subspace, i.e. r = {a — 0.5a)/{b — 
0.5(t). The system is illustrated in Fig. [T] 

The degree to which the polymer has translocated 
across the nanopore into the cavity is quantified using 
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FIG. 1. Illustration of the system studied in this work. The 
polymer translocates through a cylindrical pore of length L 
and diameter D between an ellipsoidal cavity of semi-axis 
lengths of a along z and b along x and y. The near-ellipsoidal 
space accessible to the monomer centers is shown bounded by 
the dashed blue line. 

a translocation coordinate, s, which is defined in a man¬ 
ner similar to that in a recent study of ours 4^ It is equal 
to the number of bonds that have crossed the mid-point 
of the nanopore. Typically, one bond spans this point for 
any given configuration, and this bond contributes to s 
the fraction that lies on the cavity side of the point. This 
is determined by the z coordinates of the monomers con¬ 
nected by this bond. Note that s is a continuous variable 
in the range s G [0, TV — 1]. 

III. FOKKER-PLANCK FORMALISM 


passage times is given by the probability flux at s = Sb, 

i.e. 

P(r) = J(s6,r;so,0), (3) 

where W and J are calculated using adsorbing boundary 
conditions, i.e. >V(sq, t; sq, 0) = >V(sb, t; sq, 0) = 0, and 
where the initial condition implies W = S(s — so) at t = 0. 
This represents the probability per unit time that the 
polymer reaches the boundary at s = S{,. The probability 
that the system reaches Sf, first at any time is given by 

pOO 

Pb= P{T)dT, (4) 

Jo 

and the normalized probability distribution is given by 

g{T;so) = P{T)/pb. (5) 

Finally, the associated mean first passage time is given 
by 

nOO 

(r) = / drrg(r;so). (6) 

Jo 

Similar relations can be obtained to describe transloca¬ 
tion to the other boundary. 

In this study, we use free energy functions obtained 
from MC simulations to solve Eq. m and ([2]) subject to 
absorbing boundaries at s = 0 and s = N — 1. This 
is used to calculate the mean first passage time distribu¬ 
tions in Eqs. m and ®, as well as the mean first passage 
time using Eq. and the translocation probability us¬ 
ing Eq. We use numerical integration methods, as 
described in the following section. 


We define >V(s,t;so,0) as the probability that a 
translocating polymer has coordinate s at time t given 
that it started at sq at time t = 0. In the quasi-static 
limit, yV(s,t;so,0) is governed by the Eokker-Planck 
(FP) equation. In the case where the pore friction is 
sufficiently strong the equation has the form>i 

5>V(s,t;so,0) _ 5J(s,t;so,0) 

dt ~ ds ^ ’ 

where the probability flux, so,0) is defined 


J(s,t;so,0) = -V 


( 1 dW 

\kBT ds^J ds 


( 2 ) 


where is Boltzmann’s constant and T is temperature. 
Note that the translocation rate constant 1) is indepen¬ 
dent of s in this limit. 

Consider a polymer with a translocation coordinate so 
at time t = 0, where sq lies in the domain sq G [sa,Sb]. 
For t > 0, the polymer undergoes stochastic motion gov¬ 
erned by T) and F{s) and eventually reaches the domain 
boundary at s = s^ or s = Sf,. The first passage time, 
T, is the time taken to reach s = for the first time 
without ever reaching s = Sq. The distribution of first 


IV. METHODS 

Monte Carlo simulations employing the Metropo¬ 
lis algorithm and the self-consistent histogram (SCH) 
method^ were used to calculate the free energy func¬ 
tions for the polymer-nanopore model described in Sec¬ 
tion cn The SCH method provides an efficient means to 
calculate the equilibrium probability distribution P{s), 
and thus its corresponding free energy function, F{s) = 
—fesTln'P(s). We have previously used this procedure to 
measure free energy functions in other simulation stud¬ 
ies of polymer translocatior42i2^4^ as well in a study of 
polymer segregation under cylindrical confinement.'^ 

To implement the SCH method, we carry out many 
independent simulations, each of which employs a unique 
“window potential” of a chosen functional form. The 
form of this potential is given by: 

{ oo, s < 

0, < s < (7) 

00 , s > sf^^^ 

where sf"" and are the limits that dehne the range 
of s for the i-th window. Within each “window” of s, a 
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probability distribution Pi{s) is calculated in the simu¬ 
lation. The window potential width, As = s™“ — s™™, 
is chosen to be sufficiently small that the variation in F 
does not exceed a few ^bT. Adjacent windows overlap, 
and the SCH algorithm uses the Pi{s) histograms to re¬ 
construct the unbiased distribution, V{s). The details 
of the histogram reconstruction algorithm are given in 
Ref. A description for an application to a physical 
system comparable to that studied here is presented in 
Ref. m. 

Polymer configurations were generated carrying out 
single-monomer moves using a combination of trans¬ 
lational displacements and crankshaft rotations. The 
trial moves were accepted with a probability Pacc = 
min(l, where AE is the energy difference be¬ 

tween the trial and current states. Initial polymer con¬ 
figurations were generated such that s was within the 
allowed range for a given window potential. Prior to 
data sampling, the system was equilibrated. As an il¬ 
lustration, for a. N = 121 polymer chain, the system 
was equilibrated for typically ^ 10^ MC cycles, following 
which a production run of ^ 10® MC cycles was carried 
out. During each MC cycle a move for each monomer is 
attempted once, on average. 

The windows are chosen to overlap with half of the 
adjacent window, such that = s™ 2 - The window 

width was typically As = a. Thus, a calculation for 
N = 121, where the translocation coordinate spans a 
range of s G [0,120], required separate simulations for 
239 different window potentials. For each simulation, 
individual probability histograms were constructed using 
the binning technique with 10 bins per histogram. 

The free energy functions obtained from the MC sim¬ 
ulations were used with the procedure summarized in 
Sec. lIIIl to study the translocation dynamics under the as¬ 
sumption that conformational quasi-equilibrium is main¬ 
tained during the process. The translocation probability 
W(s, t] So, 0) was determined by solving Eq. ([T|) for a cho¬ 
sen value of sq . Typically, we used sq = L/{2a), at which 
point the first monomer is just on the verge of exiting 
the pore into the cavity. The equation was solved using 
standard numerical methods with a “spatial” grid size of 
As = 0.01 and a time increment of At = 0.002'D~^. The 
distribution P(r) was calculated using Eq. ([3]), where 
a standard five-point method was used to evaluate the 
derivatives in Eq. ©■ The translocation probability, pb, 
and mean first passage time, (r), were calculated by nu¬ 
merical integration of Eqs. (IH) and (II]), respectively, using 
Simpson’s rule. 

In the results presented below, quantities of length are 
measured in units of a, energy in units of feeT, force in 
units of k^T/a and time in units of 


V. RESULTS 

We consider first the case of spherical cavities in the 
absence of a driving force or adsorption potential, i.e. 


r = 1, /d = 0 and e = 0. Figure [2] shows free energy 
functions for cavity volumes of E=150, 250 and 500, each 
for polymer lengths ranging from N=31 to A=141. As 
expected, the free energy cost of confining the polymer in 
the cavity increases as the confinement volume decreases. 
In addition, the curves all have positive curvature every¬ 
where except at near the upper and lower bounds. This 
indicates that free energy cost of inserting each monomer 
into the cavity increases as the number (and hence den¬ 
sity) of monomers inside increases. This feature is con¬ 
sistent with results from previous MC studie a^^i®^’®® and 
theoretical studie a^^i^® that explicitly account for repul¬ 
sion between monomers, in contrast to the case for ideal 
polymersj^ Another noteworthy feature is the strong de¬ 
gree of overlap between the curves for different N and the 
same V. Thus, the free energy cost of transferring one 
monomer from the outside to the inside of the cavity de¬ 
pends approximately only on the density of monomers in 
the cavity. Deviations from this trend are evident where 
F dips abruptly near s = N — 1. 



FIG. 2. Free energy functions for several different polymer 
lengths for a cavity aspect ratio of r = 1. Results for three 
different cavity volumes are shown. 

To provide a quantitative analysis of the free en¬ 
ergy functions, we employ a standard scaling theory 
approach;^ We first recall that the standard form of the 
entropic free energy barrier for translocation of a polymer 
through a narrow pore in a flat wall of negligible thick¬ 
ness is Fo/kT — (1 — A) ln[m{N — m)], where m segments 
lie on one side of the pore, and A — m lie on the other 
and where A = 0.69 is a critical exponent^ This expres¬ 
sion can easily be modified to account for the case of a 
channel in a wall of finite thickness4^ To account for the 
effect of confinement of the monomers inside the cavity, 
we consider the case where this part of the polymer is in 
the semi-dilute regime. Here, the confined portion of the 
polymer can be viewed as a collection rib blobs, each of 
size ^ ^ where v is the Flory exponent. It 

is easily shown that the confinement free energy scales 
aiM AFjkT = Ub ~ 
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m is the number of monomers in the cavity. For a finite 
length nanochannel that is spanned by an average of rip 
bonds, we note that m = s — For the pore length 

L = 1.3 used here, we estimate Up = 2 when the pore 
is filled. Thus, the total free energy, F, is expected to 
satisfy 

(F - Fp)/kT ^ (y/a3)-i-25(s _ np/2)2-25, (8) 

where we have used = |. This approximation is valid 
when np/2 < s < N — 1 — np/2. Over most of the 
range of s, rip is constant, the logarithmic term Fq is 
negligible. Only near the upper and lower bounds of s 
is the variation of Fq appreciable. This feature, as well 
as a partial emptying of the pore, accounts for the dips 
in F near s = N — Otherwise F is dominated by the 
confinement free energy of the cavity. In the case of fixed 
y, F is predicted to increase with s independent of the 
polymer length N. In addition, the variation of F with s 
has positive curvature, and F increases with decreasing 
V. These predictions are qualitatively consistent with 
the data. 

The use of the semi-dilute approximation to estimate 
free energy of confinement in the cavity is expected to be 
valid only over a restricted range of densities. In Ref. 
it was shown that the predictions are valid only up to 
packing fractions of ^ ~ 0.15, where (p = nNa^/6V. 
Beyond this value the number of monomers per blob is 
unacceptably low. A lower limit on p is imposed by the 
condition that the number of blobs, rib = be 

sufficiently large. For the polymer lengths considered in 
this work, it is difficult to find a range of s that satisfies 
both conditions simultaneously. To analyze the data, we 
use a more relaxed condition for low density and consider 
the case where rib > 3. 

Figure |3] shows the results of fits using Eq. (jS]) for 
N = 121 and various cavity volumes. The free energy 
functions have been shifted by Fo(s), which is the free 
energy function for a flat wall that was calculated ex¬ 
plicitly by simulation. We use a fitting function of the 
form F — Fq = cq -b ci(s — 1)““. The lower limit of the 
range of the fit is indicated in the plot, while the fitting 
function is deliberately extrapolated beyond the upper 
limit of the fitting range to illustrate the divergence of 
the prediction from the simulation results at high den¬ 
sity. The upper limit itself is explicitly labeled for two 
of the functions. Note that scaling predictions underes¬ 
timate the free energy in the region where p >0.15. This 
is consistent with the observation noted in Ref. that 
the confinement free energy crosses over into a concen¬ 
trated regime where the excluded volume interactions are 
screened, which leads to a higher value of a. From the fit 
to the data in the valid range, it was found that a=1.9, 
2.0 and 1.9 for E=150, 250 and 500, respectively. This 
is somewhat below the predicted value of a = 2.25. Not¬ 
ing the expected dependence of ci oc where /3=1.25, 
the ratio of ci measured for t/=500 and E=250 yields 
(3 = 0.9, while the ratio for E=500 and E=150 yields 
(3=1.1. Thus, the measured values of the exponents a 


and (3 both deviate from the predicted values. Undoubt¬ 
edly, this is arises from failing to properly satisfy the 
condition that iib >> 1- We speculate that better sat¬ 
isfying this condition would yield improved agreement. 
However, this would necessitate using polymer chains at 
least an order of magnitude larger than those considered 
here, which is not feasible for us at present. 



FIG. 3. Free energy functions for polymers of length N = 121 
and a spherical cavity of various volumes. The dashed black 
lines show fits to the data in the region for which ns > 3 and 
(j> < 0.15. The lower limit of this range is evident from the 
minimum s for the fitting curves, and the upper end of this 
range is explicitly labeled for two of the curves. For U=500, 
the upper limit extends beyond the range of the data. 

Next we consider the effects of the anisometry of the 
confining cavity. Figure S] shows free energy functions for 
chains of length N = 101 in a cavity of volume V = 500. 
Figure m^a) shows results for prolate ellipsoidal cavities 
(r > 1) and Fig. HJ^b) shows results for oblate cavities 
(r < 1). In each case, the result for spherical cavities 
is also shown, for comparison. The most notable feature 
here is the fact that deviations from spherical symmetry 
in either direction lead to an increase in the free energy. 
Note, however, that the curves for oblate cavities follow 
the opposite trend at low s (see the inset of Fig. HKb)), 
in contrast to the case for prolate cavities. The explana¬ 
tion for this difference is straightforward. As the first few 
monomers enter the cavity, the effects of confinement are 
felt mainly by the curvature of the cavity wall near the 
pore. As the cavity becomes more prolate, this local cur¬ 
vature increases, reducing the number of chain configura¬ 
tions, and F increases, accordingly. By contrast, as the 
cavity becomes increasingly oblate, the local curvature 
decreases, leading to slight initial reduction of F with 
decreasing r evident in the figure. As more monomers 
enter the cavity the polymer feels the presence of the 
cavity wall on the side opposite from the pore, and the 
trend reverses. 

Figure [5Ka) shows the variation of AF with the cav¬ 
ity aspect ratio r, where AF = F{N — 1) — F(0). Re¬ 
sults are shown for three different cavity volumes. For 
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FIG. 4. Free energy functions for a polymer of length = 101 
and a ellipsoidal cavity volume of V = 500. No driving force 
is present. Results for several different cavity aspect ratios 
are shown. The graph in (a) shows results for prolate cavities 
and (b) shows results for oblate cavities. For comparison, the 
function Fo for planar wall geometry is shown in (a). The 
inset for (b) shows a close-up of the curves at low s. 


V = 500, is a minimum for spherical cavities, as was 
noted from the results of Fig. |4l In addition, AF is ap¬ 
proximately symmetrical about r = 1 in the sense that 
AF{r) ss AF{l/r), though oblate side is slightly higher. 
For a lower cavity volume of F = 250, the trend persists, 
though with a weaker dependence of AF with r. For 
the lowest cavity volume considered, V = 150, there is 
negligible variation of AF with r. Thus, the anisometry 
of the cavity affects the confinement free energy in the 
cavity only at moderate densities, while at high densities 
no effect is observable. 

Figure EJb) shows mean hrst passage translocation 
times, (r), for ejection of the polymer from the cavity. 
The results were obtained employing the FP formalism 
described earlier, using the calculated free energy func¬ 
tions. As expected, (t) is a maximum for r=l for F=500 
and 250, and shows the same approximate symmetry 
about r=l as observed in (a). The dependence of (r) 
with r does lessen slightly as V decreases. Surprisingly, 
however, the maximum of (r) at r=l remains apprecia¬ 
ble even at F = 150, in spite of the corresponding result 



FIG. 5. (a) Free energy difference, AF = F(N — 1) — F(0), 
vs ellipsoidal cavity aspect ratio, r. Results are shown for 
A=101 for three different cavity volumes. No driving force is 
present, (b) Mean translocation time, (r), vs r. 

for AF. To understand this apparent discrepancy, let us 
compare the free energy functions for r=l and r=3 at 
each of the volumes considered. The curves are shown 
in Fig. While the variation in AF with r decreases 
with cavity volume, the free energy remains appreciably 
higher for the prolate cavity over most of the range of 
s. This corresponds to a somewhat higher average cur¬ 
vature of F for the spherical cavity. Evidently, it is this 
feature that leads to the observed difference in (r). Fig¬ 
ure [6jb) shows the full translocation time distributions 
obtained from the FP formalism, where the persistence of 
the shift toward slower translocation for aspherical cav¬ 
ities for smaller cavity volumes is clear. Note that the 
reduction of the ejection time with increasing cavity ani¬ 
sometry is consistent with the results of the Langevin 
dynamics study of Ref. [l^. 

Next, we examine the effect of the variation of AF 
with the cavity asymmetry upon changes in the polymer 
length. For a polymer at s=N — 1, there are N — 2 bonds 
inside the cavity (one is still confined in the pore). Thus, 
AF/(N — 2) is the confinement free energy per bond in¬ 
side the cavity at this point. The inset of Fig. [7] shows 
the iV-dependence of AF/{N — 2). Results are shown for 
three different cavity volumes, each with r=l and r=3. 
Several trends are evident. First, for a cavity of fixed vol¬ 
ume, the free energy per bond increases with N, and it 
decreases with F. This follows from the fact that longer 
chains and smaller volumes correspond to higher pack¬ 
ing densities inside the cavity. The other notable trend 
is that the confinement free energies for the prolate and 
spherical cavities are appreciably different at low N but 
tend to converge at high N. The convergence is stronger 
for smaller F. To clarify this trend. Fig. [7] shows the 
ratio of AF for r=3 and r=l plotted as a function of (p, 
the packing fraction in the cavity at s = N — 1. Though 
the curves for each F do not overlap, they do follow the 
same overall pattern; the ratio decreases monotonically 
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FIG. 6. (a) Free energy vs functions for a polymer of length 
N = 101. Data are shown for three different cavity volumes, 
each for spherical (r = 1) and prolate (r = 3) shapes. In all 
cases, the driving force is zero, (b) Translocation time dis¬ 
tributions calculated by numerically solving the FP equation 
and employing the free energy functions shown in (a). 


toward unity as (ft increases. The same trend was found 
for oblate (r = 1/3) cavities, as well (data not shown). 
Thus, the anisometry of the cavity has negligible effect 
on the confinement free energy at high density. 

As noted in Ref. the confinement free energy of 
a polymer under triaxial confinement scales as ^ 
semi-dilute regime, which crosses over to 
AAc ^ at higher density. In neither case does the 
expression depend on the shape of the confining cavity. 
However, these approximations likely require sufficiently 
large volumes to apply. In the former case, the num¬ 
ber of blobs is assumed to be large, implying that the 
blob dimension ^ is small relative to the confinement di¬ 
mensions. As (/ decreases, ^ increases, and the required 
condition becomes more poorly satisfied. This, in part, 
may explain the trend observed in Fig. 0 Clarification 
of these effects will require results from simulations em¬ 
ploying much longer polymer chains. 

Next, we consider the effect of a driving force acting on 
monomers inside the pore to push the polymer into the 
cavity. These calculations do not require additional sim- 



FIG. 7. Ratio of AF for r = 3 and r = 1 vs maximum packing 
fraction, (j). Data are shown for polymer lengths in the range 
A=31-141 and cavity volumes of V =150, 250 and 500, each 
with aspect ratios of r = 1 and 3. The inset shows the scaled 
free energy barrier heights, AF/{N — 2), vs polymer length, 
N. Solid symbols correspond to r = 1 and open symbols 
correspond to r = 3. 

ulations. To an excellent approximation, the free energy 
function for finite /d is given by 

F{s;U)=F{s-0)- ULs (9) 

Thus, previous results for /d = 0 can be easily modified 
to yield free energy functions for finite /d. A few simu¬ 
lations with a finite driving force were carried out, and 
the validity of this expression was confirmed (data not 
shown). 

To study the effect of varying V and r on F{s) and the 
translocation rate, it is helpful to introduce a modified 
driving force, defined /^ = /d - AF(/d = 0)/[L{N - 1)]. 
Note that /^ = 0 corresponds to the case where F{Q) = 
F{N — 1) for arbitrary /d and V. The shift is designed 
to examine the combined effect of the driving force and 
the curvature of F{s) in the absence of the entropic bias 
toward polymer ejection. 

Figure [8] shows results for polymers of length N = 101 
in a cavity of volume F=150 and F=500. Data are 
shown for several different driving forces, each for the 
case of a spherical cavity and a prolate (^=3) cavity. 
The free energy curves in Fig. [SKa) are also labeled with 
f^. The curves in Fig. HKb) for F=500 correspond to 
specific values of //. For F=500, note that r=l and 
r=3 curves with the same correspond to different /d, 
since AF(/d = 0) differ for different r, as was evident 
in Fig. EKa). The main effect of the driving force is to 
remove the free energy penalty for inserting the poly¬ 
mer into the cavity. This point is first reached when 
/^=0 (i.e. /d = AF(/d = 0)/[L{N - 1)].) Beyond that 
point F is lower for the polymer completely inside than 
where it is completely outside. As shown for F=150 in 
Fig. [i(a), increasing /d generally increases rate that F 









decreases with s. However, the magnitude of dF/ds de¬ 
creases with s, which is a consequence of the curvature 
of F{s) for /d=0. In addition, for moderate values of /d, 
the function exhibits a local minimum when the polymer 
is mostly inside the cavity. For example, for /d=1.3, the 
minimum lies at Smin ~80, as is clear from the inset of the 
figure. The magnitude of the resulting barrier between 
Smin and s = N — \ decreases and eventually vanishes 
about some critical value of /d. The presence of such a 
minimum has been noted in previous studies of translo¬ 
cation into spherical cavities with a driving force arising 
from attraction of monomers to the cavity walls and 
effective forces arising from translocation between cavi¬ 
ties of different sizes 



FIG. 8. Free energy functions for a polymer of length A'^=101 
and a cavity of volume F=150 in (a) and y=500 in (b). Re¬ 
sults are shown for several /d and/or where the shifted 
driving force is defined in the text. Solid curves corre¬ 
spond to spherical cavities (i.e. r=l) and dashed curves are 
shown for prolate cavities with r=3. The inset in (a) shows a 
close-up of the curves for /d=1.3. 

As noted earlier, AF {= F{N — 1) — F(0)) is approxi¬ 
mately invariant with respect to r for N=101 and F=150 
in the case where /d=0. From Eq. Q it follows that this 
is also true for arbitrary /d. However, over most of the 
range of s, there is an appreciable dependence of F{s) 
on r. This is evident in Fig. [SJa), which illustrates the 
difference in the curves for the cases of r=l and r=3. 


Unlike the case illustrated in Fig. |6]for /d=0, where the 
anisometry of the cavity has a small effect on the rate for 
polymer ejection, small changes in r can have a signifi¬ 
cant effect on the rates of polymer insertion. 



FIG. 9. Mean translocation time (r) vs shifted driving force, 
/d, where /d = /d — AF{fd = 0)/{L{N — 1)). Results are 
shown for a polymer of length A=101 for cavity volumes of 
U=150 and U=500, each for three different cavity aspect ra¬ 
tios. The inset shows the same results plotted as a function 
of the driving force, fd- 

The effects of the free energy functions on the polymer 
insertion translocation times are illustrated in Figure IHl 
which shows results for translocation times calculated us¬ 
ing the FP method for Af=101 and volumes of V =150 and 
U=500, each for cavity anisometries of r=l/3, 1, and 3. 
The inset shows (r) vs fd- As expected, the translo¬ 
cation times decrease with increasing fd- In addition, 
the translocation rate is faster when the cavity volume is 
larger. This follows from the fact that the entropic force 
acting against fd is greater for smaller volumes, as was 
noted in Fig. Ha). The main part of the figure shows (t) 
vs /d. At large /j, the translocation times all converge to 
the same curve. In this regime, effects of the curvature 
of F{s) become negligible, and the translocation rate is 
governed principally by dF/ds- For any given value of 
dF/ds ~ fd independent of V, and so (r) will be equal. 
At lower values of Z^, the results for the two volumes 
diverge. The curves for U=150 begin to increase rapidly 
near RiO.5 and those for I/=500 diverge near ~0.1. 
In each case, the rapid increase in (r) is due to the in¬ 
creasing influence of the curvature of F{s) as the driving 
force decreases. The curvature gives rise to the free en¬ 
ergy minimum evident in Fig. |S]and the associated free 
energy barrier. As was clear in Fig.H^')) the curvature of 
the function is greater when the cavity volume is lower. 
Thus, the onset of the rapid increase in (r) occurs first 
for U=150. For this volume, the free energy minimum 
first appears around fd ~ 1.5, and deepens rapidly as fd 
increases. As seen in Fig.H this corresponds to ~ 0.7. 



















9 


By contrast, the free energy minimum for V =500 appears 
at /j ssO.l. Both of these results are consistent with the 
trends in the figure. 

Finally, we consider the effects of cavity anisometry 
on the translocation times. At high /j, it is clear that 
varying r has negligible effect on (r). In this limit, 
F{s) ~ —fdLs, which is independent of r, and so the 
invariance of (t) to r is expected. For F=150, a dif¬ 
ferent trend emerges at lower /d, where the free energy 
minimum emerges. In this regime, (t) becomes greater 
for spherical cavities than for either prolate or oblate 
cavities. For example, for /d = 1.15, (r) is 1.6 times 
greater for r=l than for r=3. This effect is not present 
for l/=500. Thus, shape anisometry of the cavity leads 
to faster polymer insertion at low driving force and high 
packing fraction. This can be understood from inspec¬ 
tion of Fig. |Sl^a), where we observe a deeper free energy 
minimum for r=l than for r=3. This trend persists for 
functions with lower /d than those values shown in figure. 
In the case of F=500, there appears to be negligible ef¬ 
fect of the anisometry on (r). From the arguments above, 
there is expected to be some effect in the regime where 
a free energy minimum is present. From Fig. |Sl[b), this 
occurs for ^0.1 However, we note that the difference 
in the depths of the minima for r=l and r=3 is < kT. By 
contrast, the difference is considerably greater for F=150 
at some positive values of /^. Consequently, the effects 
of the cavity anisometry are expected to be more appre¬ 
ciable for smaller F, in accord with the results. 

The translocation time predictions can be compared 
with the results of the Langevin dynamics simulation 
study of Ali et al^ They studied ejection and insertion of 
flexible and semi-flexible polymers of length N = 100 for 
a prolate ellipsoidal cavity of comparable dimensions to 
the V = 150 cavity considered here, as well as for a spher¬ 
ical cavity of the same volume. In the case of ejection of 
a flexible polymer, they found faster translocation for el¬ 
lipsoidal cavities, in agreement with the results presented 
earlier. However, for insertion of a flexible polymer, they 
found that translocation was faster for the spherical cav¬ 
ity. Those results are in disagreement with our predic¬ 
tions. Most likely, this is due to out-of-equilibrium con¬ 
formational behaviour in the dynamics simulations. Such 
effects are not accounted for in the FP approach, which 
assumes quasistatic dynamics. The predictions are also 
in disagreement with the 2-D Langevin dynamics study 
in Ref. l3lL where insertion proceeded most rapidly in the 
case of a circular cavity, independent of packing fraction. 
In this case, it is less clear whether the discrepancy arises 
from nonequilibrium dynamics or from a different under¬ 
lying behaviour of the free energy for 2-D. Regardless, 
these results demonstrate the need to exercise caution 
in using free energy arguments to interpret translocation 
simulation results. The assumptions for the valid use of 
the FP approach are only satisfied for sufficiently high 
pore friction, as we noted in Refs. HI and However, 
such predictions are useful for identifying some of the 
qualitative effects of nonequilibrium dynamics. 


The final case we consider in this study is transloca¬ 
tion into an adsorbing cavity. Here, monomers have an 
energy of —e if their centers lie within a distance a (=1) 
from the elliptical surface of the cavity. Figure [TUI shows 
free energy functions for various values of e for a poly¬ 
mer of length A^=101. Results are shown in (a) for a 
cavity volume of F = 150 and in (b) for F=500, each for 
cavity anisometries of r=l and 3. As expected, increas¬ 
ing the attraction to the cavity wall decreases AF, and 
thus removes the entropic cost of insertion of the polymer 
into the cavity. Generally, the free energy function has 
positive curvature. Beyond some value of e, we see that 
AF < 0, i.e. it is more favourable for the polymer to 
lie completely inside rather than completely outside the 
cavity. However, there is an intermediate range of e for 
which there is a local free energy minimum and thus a free 
energy barrier that must be overcome for complete inser¬ 
tion. The depth of the free energy minimum is greater 
in the cavity with the smaller volume. These trends are 
qualitatively comparable to those for the curves shown 
in Fig. [5] for the case of a driving force in the nanopore. 
In addition, the curves a re q ualitatively similar to those 
calculated in Refs. andl27l for adsorbing spherical cavi¬ 
ties. Note that the details of the model employed in those 
studies differed from the present model, and thus a direct 
quantitative comparison of the results is not possible. 

For each e value, the free energy curve for the pro¬ 
late cavity lies above that for the spherical cavity. The 
results for oblate cavities with r=l/3 were comparable 
to those of r=3 (data not shown). Thus, increasing the 
anisometry of the cavity reduces the effectiveness of the 
attraction to the cavity wall to drive the polymer into 
the cavity. This general result is also clear in the inset 
of Fig. fTOT al. which shows AF vs e. Comparing with 
Fig. [51 we note that the cavity anisometry has the same 
effect on the free energy for the case of attractive cavity 
walls as it did for the case of a driving force acting in the 
pore. However, the differences between the free energies 
for spherical and anisometric cavities are greater in the 
present case, and increase with increasing e. Comparing 
the results in Fig. m a) and (b), we find that the abso¬ 
lute degree to which the the curves are shifted between 
r=l and r=3 is not significantly affected by the cavity 
volume. 

Translocation time distributions calculated using the 
FP formalism and the free energy functions in Fig. fTHT bl 
are illustrated for A=101 and F=500 in Fig.[Tl] Results 
are shown for values of e that are sufficiently large for an 
appreciable probability that the polymer is driven into 
the cavity. As expected from the free energy functions, 
increasing the cavity anisometry significantly slows down 
the rate of insertion. By contrast, cavity anisotropy has 
a negligible effect on the polymer insertion rate for cav¬ 
ities of the same volume in the case where insertion is 
driven by a force acting in the nanopore, as was evident 
in Fig. ini The strong effect observed here apparently 
arises from the fact that changing the shape of the cav¬ 
ity distorts the polymer in a way that reduces the average 
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FIG. 10. Free energy functions for an adsorbing cavity for 
various adsorption strengths for a polymer of length A^=101. 
Results in (a) are for F^ISO and (b) are for 1/=500. The 
solid curves are for spherical (r=l) cavities, and the dashed 
curves are for prolate {r=3) cavities. The inset in (a) shows 
the free energy difference AF = F{N — l)-F’(O) vs e. 


number of monomer-surface contacts in the cavity for the 
accessible conformations. For the model with the driving 
force in the nanopore, this kind of effect does not affect 
the energetic contribution to the free energy, unlike the 
case where surface interactions are present. 


VI. CONCLUSIONS 


In this study, we have used simulation and theoreti¬ 
cal methods to investigate the translocation of a flexible 
hard-sphere polymer into and out of an ellipsoidal cav¬ 
ity. Monte Carlo simulations were employed to calculate 
translocation free energy functions, and these functions 
were used together with the Fokker-Planck equation to 
predict the translocation times. We considered the case 
of ejection from the cavity for an athermal system, as 
well as the case of polymer insertion driven either by a 
force acting in the pore or by monomer attraction to the 
walls of the cavity. We studied the effects on the free 
energy and predicted dynamics of varying all of the rel- 



FIG. 11. Translocation time distributions for systems with 
an attractive cavity wall, for A^=101 and F=500. Results for 
two values of e are shown, each for cavity aspect ratios of r=l 
and r—3. 


evant system parameters, with special attention to the 
effect of the cavity anisometry. 

The free energy functions calculated for the ather¬ 
mal system with spherical cavities are consistent with 
those calculated previously by analytical methodsi^ and 
by simulationi^L^ We find that the variation of F with 
packing fraction is approximately consistent with a scal¬ 
ing theory prediction that treats the polymer to be in the 
semi-dilute regime inside the cavity. Deviations from the 
predicted scaling exponents likely arise from finite-size 
effects. Increasing the anisometry of the cavity at fixed 
volume increases the free energy cost of insertion. The 
effect is approximately symmetric with respect to pro¬ 
late or oblate distortions and diminishes with increasing 
cavity packing fraction. This increase in the free energy 
leads to a prediction of faster polymer ejection for el¬ 
lipsoidal cavities compared to spherical cavities. This is 
consistent with the ejection rates measured in a previous 
Langevin dynamics simulation studyi^ Application of a 
driving force on monomers in the pore leads to a local 
free energy minimum for an intermediate range of force 
magnitudes, and the corresponding free energy barrier 
slows the insertion. Increasing the shape anisometry of 
the cavity causes faster polymer insertion at high packing 
fractions and moderate driving force, but has negligible 
effect for larger cavities or high driving force. This re¬ 
sult is inconsistent with Ref. [^, most likely because of 
out-of-equilibrium conformational behaviour that is not 
accounted for in the FP approach. Finally, the translo¬ 
cation rate for insertion driven by monomer attraction to 
the cavity walls is much more sensitive to cavity anisom¬ 
etry than for the case of insertion driven a force in the 
pore. 

It is intriguing that the predicted qualitative effects of 
the cavity shape on translocation rate are in agreement 
with observations from dynamics simulations for ejection 
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but not insertion. As noted, the latter discrepancy is dne 
to nonequilibrium behaviour, most likely associated with 
the encapsulated part of the polymer, which provides 
the dominant contribution to the free energy function. 
The internal dynamics of this portion likely becomes very 
sluggish at high densities. Increasing the cavity anisom- 
etry at fixed volume narrows the cavity in at least one 
dimension, leading perhaps to slower conformational re¬ 
arrangement and an increase in nonequilibrium effects. 
However, it seems likely that the ejection process would 
also exhibit this behaviour as well for the same simula¬ 
tion model. Evidently, there is a hysteresis effect at play 
here, which was also noted in Ref. [2^ In future work, 
we will examine this issue further and quantify the devi¬ 
ation from equilibrium behaviour by carrying dynamics 
simulations for variable pore friction. At high pore fric¬ 
tion, the polymer ejection dynamics are perfectly con¬ 
sistent with FP predictions;^ and it is likely to be the 
case for insertion as well. This work should contribute 
to understanding the limits of the validity of employ¬ 
ing free energy arguments to understand DNA ejection 
and packaging in viruses, as well as other systems, where 
nonequilibrium effects are appreciable. 
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